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Abstract 

We theoretically investigate the spatial dependence of the order parameter of 
the two-dimensional semi-infinite Ising model with a free surface at or above 
the bulk critical temperature. Special attention is paid to the influence of 
a surface magnetic field h\ and the crossover between the fixed points at 
h\ = and hi = oo. The sharp increase of the magnetization m(z) close 
to the boundary generated by a small hi, which was found by the present 
authors in the three-dimensional model, is also seen in two dimensions. There, 
however, the universal short-distance power law is modified by a logarithm, 
m(z) ~ z re lnz, where k, the difference between the scaling dimensions of 
hi and the bulk magnetization, has the exact value 3/8. By means of a 
phenomenological scaling analysis, the short-distance behavior can be related 
to the logarithmic dependence of the surface magnetization mi on h\. Our 
results, which are corroborated by Monte Carlo simulations, provide a deeper 
understanding of the existing exact results concerning the local magnetization 
and relate the short-distance phenomena in two dimensions to those in higher 
dimensionality. 
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I. INTRODUCTION 



In systems at or close to a critical point the asymptotic power laws governing the behavior 
of thermodynamic quantities are modified in the vicinity of surfaces or other imhomogenities 
HIIJ]. The characteristic distance within which changes occur is given by the bulk correla- 
tion length £. In general each bulk universality class splits up in several surface universality 
classes, depending upon whether the tendency to order in the surface is smaller or larger 
than or the same as in the bulk. In the case of the two-dimensional (2-d) Ising model in 
a semi-infinite geometry there exist two surface universality classes. Since the boundary is 
one-dimensional and, thus cannot become critical itself, the surface generally reduces the 
tendency to order. More precisely, for any positive starting value the exchange coupling 
between surface spins J\ is driven to zero by successive renormalization-group transforma- 
tions. The relevant field pertaining to the surface is the magnetic field hi, which acts on 
surface spins only and which for instance may take into account the influence of an adjacent 
noncritical medium. The two universality classes then are labelled by hi = 0, the "ordinary" 
transition, and hi = oo, the "normal" transition, where the former is a unstable fixed point 
and the latter is a stable fixed point of the renormalization-group flow PL0]. 

In this work we are mainly concerned with order-parameter profiles for finite hi in the 
crossover region between the fixed points. Nonetheless, let us first recapitulate the situation 
at the fixed points. For the sake of simplicity, in the Introduction our considerations remain 
restricted to bulk criticality T = T c , but the extension to the critical region is straightforward 
and will be done below. At the critical point and hi = 0, the order-parameter (or magneti- 
zation) profile m(z) is zero for any distance z > from the surface. In the other extreme, 
hi = oo, it is well known that at macroscopic distances from the surface the magnetization 
decays as ~ z~ X4 ', where = /3/u is the scaling dimension of the bulk order-parameter 
field, with the exact value 1/8 in the 2-d Ising model. 

What do we expect, when hi takes some intermediate value, i.e., in the crossover region 
between the fixed points? Now, hi will certainly generate a surface magnetization mi. As 
far as the profile m(z) is concerned, a first guess would perhaps be that the magnetization 
should decay from that value as z increases away from the surface. This guess is supported by 
mean-field theory, where one can calculate m(z) and indeed finds a monotonously decreasing 
function of z 0j|[J. 

In Ref. || the present authors have shown that, contrary to the naive (mean-field) expec- 
tation, fluctuations may cause the order parameter to steeply increase to values m(z) ^> mi 
in a surface-near regime. The range within which this growth occurs (at bulk criticality) is 
determined by hi, the characteristic length scale being l\ ~ h' x where yi — Ai/u is the 
scaling dimension of hi 0. As further demonstrated by the present authors in Ref. ||, the 
growth of order in the near-surface regime z <C li is described by a universal power law 

m ~ z K with k = yi — x<f, , (1) 

i.e., the growth exponent k is governed by the difference between the scaling dimensions 
yi and x^. For z ~ li the profile has a maximum and farther away from the surface, at 
distances much larger than li, the magnetization decays as z~ x '*' . 
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Largely analogous results — monotonous behavior at the fixed points and profiles with one 
extremum in the crossover regime — where previously found by Mikheev and Fisher [§ for 
the energy density of the 2-d Ising model. The authors also suggested to calculate the order 
parameter in the crossover region. Below we focus our attention exactly on this problem. 
The questions posed in this work are the following: 

• Does the short- distance growth of m(z), found for instance in the three-dimensional 
Ising model, also occur in two dimensions? 

• Does the simple power law (p]) quantitatively describe the magnetization profile in 
d = 2, or are modifications to be expected? 

• Is the scenario for the crossover between "ordinary" and "normal" transition, developed 
for the three-dimensional Ising model, also valid in d = 2? 

As we will demonstrate below, the answer to the first and third question is "yes", but the 
simple power law (JJ) is modified by a logarithm. 

The rest of this paper is organized as follows: In Sec.|TI| the theoretical framework is 
expounded. We first summarize the results of Ref. [[5] and then generalize the scaling 
analysis by taking into account the available exact results on the dependence of mi on hi 
and on the magnetization profile. In Sec.|T|, in order to corroborate our analytical findings, 
we present Monte Carlo (MC) data for m(z) and compare with exact results for the order- 
parameter profile. In two Appendices exact literature results on the dependence of the 
surface magnetization on hi and on the order-parameter profile are briefly reviewed. 



II. THEORY 

We consider the semi-infinite Ising system with a free boundary on a plane square lattice. 
The exchange coupling between neighboring spins is J. A surface magnetic field Hi is 
imposed on the boundary spins and bulk magnetic fields are set to zero such that the 
Hamiltonian of the model reads 

W = -J E s i s i - H i E s * ' ( 2 ) 

<ij>£V iedV 

where dV and V stand for the boundary and the whole system (including the boundary), 
respectively. As usual, we work with the dimensionless variables 

K = J/k B T and h x = H x jk B T . (3) 

The bulk critical point corresponds to K c = | ln(l + y/2). 
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A. Scaling analysis for Ising system in 2 < d < 4 



In the critical regime, where |r| = \(T — T c )/T c \ <C 1, thermodynamic quantities are 
described by homogeneous functions of the scaling fields. As a consequence, the behavior of 
the local magnetization under rescaling of distances should be described by 

m{z, r, h x ) ~ b~ x * m(z& -1 , rb 1/u , h x b yi ), (4) 

where = (3/v and yi = Aj/i/ is the scaling dimension of /ii 0. In terms of other surface 
exponents we have yi = (d — r}\\)/2 In Eq. @ it was further assumed that the distance 
z from the surface is much larger than the lattice spacing or any other microscopic length 
scale. One is interested in the behavior at macroscopic scales, and, for the present, z may 
be considered as a continuous variable ranging from zero to infinity. 

Removing the arbitrary rescaling parameter b in Eq. by setting it ~ z, one obtains 
the scaling form of the magnetization 

m(z, t, ~ z"^ M(z/£,z/h) , (5a) 

where, as already stated, 

h ~ h^ 1/yi (5b) 

is the length scale set by the surface field. The second length scale pertinent to the semi- 
infinite system and occurring in ( |5a| ) is the bulk correlation length £ = t~ u . Regarding 
the interpretation of MC data, which are normally obtained from finite lattices, one has to 
take into account a third length scale, the characteristic dimension L of the system, and 
a finite-size scaling analysis has to be performed. The latter will be briefly described in 



Sec. IIC 



Going back to the semi- infinite case and setting r = 0, the only remaining length scale 
is l\, and the order-parameter profile can be written in the critical-point scaling form 

m{zM)~z~ x *M c {z/h). (6) 

As said above, for z — ► oo the magnetization decays as ~ z~ x ^ and, thus, M. c {£) should 
approach a constant for C, — > oo. In order to work out the short- distance behavior of the 
scaling function Ai c (Q, we demand that m(z) ~ mi as z — > 0. This means that in general, in 
terms of macroscopic quantities, the boundary value of m(z) is not mi. If the ^-dependence 
of m(z) is described by a power law, it cannot approach any value different from zero or 
infinity as z goes to zero. However, the somewhat weaker relation symbolized by "~" should 
hold, stating that the respective quantity asymptotically (up to constants) "behaves as" or 
"is proportional to". This is in accord with and actually motivated by the field-theoretic 
short-distance expansion [|,@, where operators near a boundary are represented in terms of 
boundary operators multiplied by c-number functions. 

In the case of the three-dimensional Ising model the foregoing discussion leads to the 
conclusion that m(z) ~ hi because the "ordinary" surface — the universality class to which 
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also a free surface belongs — is paramagnetic and responds linearly to a small magnetic field 
||. The consequence for the scaling function in (H) is that M c (() ~ C yi > an d, inserting this 
in (|^), we obtain that the exponent governing the short-distance behavior of m(z) is given 
by the difference between y 1 and x$ (as already stated in Eq. ([I])). Using the scaling relation 
V± = (v + ? 7||)/2 @ among anomalous dimensions, one can reexpress the exponent k as j/J 

k = 1 - V± • (7) 

In the mean field approximation the value of k is zero, and one really has m(z — > 0) = mi. 
However, a positive value is obtained when fluctuations are taken into account below the 
upper critical dimensionality d*. For instance, the result for the n- vector model with n = 1 
(belonging to the Ising universality class) in the framework of the e-expansion is k — e/6 || . 
Thus, the magnetization indeed grows as z increases away from the surface. 

For the 2-d Ising model the exponent rj± is known exactly [[J and one obtains k = 3/8. 
However, as will be discussed in Sec. [11 BL the pure power law found in Ref. || for d — 3 is 



modified in two dimensions by a logarithmic term, and the exponent 3/8 cannot directly be 
seen in the profile. 

The above phenomenological analysis is straightforwardly extended to the case r > 0. 
In d > 2, we may assume that the behavior near the surface for z << £ is unchanged 
compared to ([[]), and, thus, the increasing profiles are also expected slightly above the 
critical temperature. The behavior farther away from the surface depends on the ratio h/!;. 
In the case of l\ > £ a crossover to an exponential decay will take place for z ~ £ and the 
regime of nonlinear decay does not occur. For l\ < £ a crossover to the power-law decay 
~ z~Pl v takes place and finally at z ~ £ the exponential behavior sets in. 

Below the critical temperature, the short- distance behavior of the order parameter is 
also described by a power law, this time governed by a different exponent, however. The 
essential point is that below T c the surface orders even for hi = 0. Hence, in the scaling 
analysis the scaling dimension of h\ has to be replaced by the scaling dimension of mi, the 
conjugate density to hi, given by x\ = (3\/u |fj]. The exponent that describes the increase 
of the profile is thus X\ — [|l(J , a number that even in mean-field theory is different from 



zero (= 1/2). 

Phenomena to some extent analogous to the ones discussed above were reported for the 



crossover between special and normal transition [pj] , ^ . Also near the special transition 
the surface field h\ gives rise to a length scale. However, the respective exponent, the 
analogy to n in (|I|), is negative, and, thus, one finds a profile that monotonously decays for 
all (macroscopic) z, with different power laws in the short-distance and the long-distance 
regime and a crossover at distances comparable to the length scale set by h±. However, 
non-monotonous behavior in the crossover region is a common feature in the case of the 
energy density in d = 2 and as well as in higher dimensionality |13|] . 

The spatial variation of the magnetization discussed so far strongly resembles the time 
dependence of the magnetization in relaxational processes at the critical point. If a system 
with nonconserved order parameter (model A) is quenched from a high-temperature initial 
state to the critical point, with a small initial magnetization m' 1 ', the order parameter 



behaves as m ~ t e , where the short-time exponent 9 is governed by the difference 
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between the scaling dimensions of initial and equilibrium magnetization divided by the 
dynamic (equilibrium) exponent [15] . Like the exponent k in ([!]), the exponent 9 vanishes 
in MF theory, but becomes positive below d*. 



B. Scaling analysis in d = 2 



During the years, initiated by the work of McCoy and Wu p6|J17| , the 2-d Ising model 
with a surface magnetic field received a great deal of attention, because many aspects can be 
treated exactly and it is a simple special version of the 2-d Ising model in a inhomogeneous 
bulk field, a problem to which an exact solution would be highly desirable. Some of these 
exact results, namely those considering the vicinity of the critical point |18"l , |T9l , will be used 
in the following as a guiding line for our phenomenological scaling analysis and to compare 
numerical data with. 

The dependence of mi on the surface magnetic field (bulk field h = 0) in the 2-d Ising 
model was calculated exactly by Au-Yang and Fisher JT8] in a n x oo (strip) geometry. The 



limit n — > oo, yielding results for the semi- infinite geometry, was also considered. Whereas 
above two dimensions at the ordinary transition the surface, in a sense, is paramagnetic, i.e., 
the response of mi to a small hi is linear, in two dimensions the function mi(/ti) has a more 
complicated form. As summarized in Appendix A (see Eq. (|A2[ )), there is a logarithmic 
correction to the linear term in d = 2; for hi — > the surface magnetization behaves as 
~ hi In hi. 

Further, as shown by Bariev [Tj| and summarized in Appendix B, the length scale li 
determined by hi behaves as ~ [tanh(/ii)]~ 2 . For small h\, where the scaling analysis is 
expected to be correct, this is consistent with (EB) as yi = 1/2 in the 2-d Ising model. Thus 
the characteristic length scale that enters the scaling analysis depends in the same way upon 
hi as in higher dimensions. 

The foregoing discussion allows us to generalize our scaling analysis, especially Eqs. @ 
and the near-surface law ([I]), such that the speacial features of the 2-d Ising model are taken 
into account. Again, the only available length scale at r = is h, and the magnetization 
can be represented in the form given in Eq. For z — > oo we expect that m ~ z -1 / 8 and, 
thus, Ai c (() should approach a constant for £ — > oo. In order to find the short- distance 
behavior we assume again that m(z) ~ mi as z —>■ 0. Taking into account the logarithmic 
correction mentioned above and discussed in more detail in Appendix A (see Eq. ( |A2| ), we 
find that Ai c (() ~ ^ 1 ' 2 ln^ for £—►(). Hence, for the short- distance behavior of m(z) in the 
semi-infinite system we obtain 

m(z,hi)~h 1 z K ki(h 1 z V1 ), (8) 

where the exact values of the exponents are k = 1 — r]± = 3/8 and yi = 1/2. Thus, for 
z < li the magnetization m(z) for a given value of hi behaves as ~ z 3 ^ 8 hiz. 

The result (§) should hold for any value of the exchange coupling Ji in the surface. 
In our MC analyses to be presented below we implemented free boundary conditions with 
J = Ji, but we expect (§) to hold for any value of Ji with possible Ji-dependent nonuniversal 
constants leaving the qualitative behavior of the profiles unchanged. 
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Eq. (H) is the main analytic result of this work. As discussed in the following, it is consis- 
tent with Bariev's exact solution Jl9| (see Appendix B) and with MC data for the profile. It 
tells us that the short- distance power law behavior is modified by a logarithmic term. This 
logarithm can be traced back to the logarithmic singularity of the surface susceptibility [|TJ, 
which, in turn, causes the logarithmic dependence of mi on hi, and eventually leaves its 
fingerprint also on the near-surface behavior of the magnetization. The result (|8|) provides 
a thorough understanding of the near-surface behavior of the order parameter and allows 
to relate special features of the two-dimensional system, which were (as we will discuss in 
more detail below) previously known from the exact analyses [|Tj|] , to the somewhat simpler 
short- distance law in higher dimensions. 



C. Finite Size scaling 



In order to assess the finite size effects to be expected in the MC simulations, we have to 
take into account the finite-size length scale L, which is proportional to the linear extension 
N of the lattice (compare Sec. [Ill A] below). The generalization of (Q) reads [SC 



m(z, r, hi, L) ~ b~ x * mizb- 1 , rb 1/u , hi b yi , Lb' 1 ) , (9) 

and proceeding as before, we obtain as the generalization of fl5a|) to a system of finite size: 

m(z, t, h, L) ~ z-*+M{z/t z/h, z/L) . (10) 

Thus even at T c there are two pertinent length scales, on the one hand L (imposed by the 
geometry that limits the wavelength of fluctuations) and on the other hand li (the scale set 
by hi). 

It is well known that for large z > L we have to expect an exponential decay of m(z) on 
the scale L. In the opposite limit, when z is smaller than both L and li, we expect the short- 
distance behavior @ to occur. However, for d — 2 it can be concluded from the finite-size 
result ( |A4j ) that there will be an L-dependent amplitude, a prefactor to the function given in 
(§). But otherwise the logarithmically modified power law ([§]) should occur. Farther away 
from the surface, the form of the profile depends on the ratio between li and L. For z smaller 
than both L and li, the behavior is described by (|8]). In the case of li > L a crossover to 
an exponential decay will take place for z ~ L. In the opposite case, a crossover to the 
power-law decay ~ z~ x ' t ' takes place, followed by the crossover to the exponential behavior 
at z ~ L. Thus, qualitatively, the discussion is completely analogous to the one in Sec. |ll A| , 
where the behavior at finite £ was described. 



III. MONTE CARLO SIMULATION 



A. Method 



The results of the scaling analysis, especially the short-distance law (||), were checked by 
MC simulations. To this end, we calculated order-parameter profiles for the 2-d Ising model 
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with uniform exchange coupling J. The geometry of our systems was that of a rectangular 
(square) lattice with two free boundaries (opposite to each other) and the other boundaries 
periodically coupled, such that the effective geometry was that of a cylinder of finite length. 
The linear dimension perpendicular to the free surfaces was taken to be four times larger 



than the lateral extension in order to keep corrections due to the second surface small pi 
Hence, when we talk about a lattice of size N in the following, we refer to a rectangular 
N x AN system. 

In order to generate an equilibrium sample of spin configurations, we used the Swendson- 
Wang algorithm JEJ. It effectively avoids critical slowing down by generating new spin 
configurations via clusters of bonds, whereby the law of detailed balance is obeyed. For 
a given spin configuration, a bond between two neighboring spins of equal sign exits with 
probability 1 — e~ 2K . There are no bonds between opposite spins. Then, clusters are defined 
as any connected configuration of bonds. Also isolated spins define a cluster, such that 
eventually each spin belongs to one of the clusters. After having identified the clusters, the 
new configuration is generated by assigning to each cluster of spins a new orientation, with 
equal probability for each spin value as long as the cluster does not extend to a surface. 

In order to take into account hi, we introduced, in the same way as suggested by Wang 
for taking into account bulk fields , two "ghost" layers of spins next to each surface that 
couple to the surface spins with coupling strength hi and that all point in the direction of 
hi. If at least one bond between a surface and a ghost spin exists the cluster has to keep 
its old spin when the system is updated. This preserves detailed balance. In the practical 
calculation this rule was realized by a modified (reduced) spin-flip probability 

p(k) = 1 - - exp(-2 hi k) (11) 

for clusters pointing in the direction of hi (and 1/2 for clusters pointing in opposite direc- 
tion), k being the number of surface spins contained in the cluster. 

In order to obtain an equilibrium distribution of configurations we discharged several 
hundred (depending on system size) configurations after the start of the simulation. To 
keep memory consume low, we used multispin-coding techniques, i.e., groups of 64 spins 
were coded in one long integer. 



B. Comparison with exact results 

A crucial test for the MC program is the comparison with known exact results. On the 
other hand, if both MC data an exact results agree, the former can be regarded also as a 
confirmation of the exact results. Having calculated order-parameter profiles for different 
values of hi, in particular the magnetization at the boundaries (the ends of the cylinder) 
can be obtained. In Fig. 1 we show the results for mi(hi). The squares respresent our MC 
data for a system of size N = 512. Also depicted is the exact result for the semi-infinite 
system taken from Ref. (see also Appendix |A|). It is clear that the MC values approach 
the exact curve for large hi. Below hi ~ 0.03, the MC results show a linear dependence 
on hi, significantly deviating from the exact curve. The linear behavior can be regarded as 
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a finite-size effect, and it is qualitatively consistent with the result (|A4|) of Au-Yang and 
Fisher [18]. The latter was derived in a strip of finite width with infinite lateral extension, 
however, so that a quantitative comparison with our data is not possible. 




0.12 



Fig. 1: Monte Carlo results for mi as a function of h\ for N = 512 (represented by full squares) 
compared to the exact result of Ref. Jl8| (see Appendix A). The statistical errors of the Monte 
Carlo data are about the same size as the symbols. A detailed discussion and comparison of the 
data is in the text. 



Next we compare with the exact solution obtained by Bariev |TJ| for the order-parameter 
profile. The explicit result ( P5| ) holds in the limit h\ — > 0. This limit is hard to access in 
the MC simulation since the signal m(z) becomes small and is eventually lost in the noise. 
To obtain a concrete result to compare with, we have calculated the profile numerically 
from (|B4|) terminating the series in ( B4a|) after the third term. It turned out that the series 
converges rather quickly as long as the distance from the surface is not too small. Only 
very close to the surface higher orders need to be taken into account. Concretely, we took 
K = 0.9997^ (i.e. r ~ 0.001) and h x = 0.01. Then, employing (BJ) and O one obtains 
£ = 567 and li = 2069, respectively. The result of the numerical evaluation of (|B4] ) is shown 
in Fig. 2 (dashed curve), where m(z) versus the distance z = n — 1 is plotted. 

Then, with exactly the same parameters, the MC profiles were calculated, with the size 
N varying between 128 and 1024. The results are depicted in double-logarithmic represen- 
tation in Fig. 2 (solid curves). It is obvious that the MC data approach the (approximated, 
in principle) exact profile with increasing lattice size. Most importantly, both results show 
the short- distance behavior anticipated from the scaling analysis above and expressed in (^j). 
This is demonstrated by the dotted line, which depicts the function 0.016 n 3 / 8 (6.15 — Inn) 
(the constants were fitted). As expected, it only describes the profile for short distances 
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from the boundary and becomes wrong for large distances, completely analogous to the 
asymptotic form (A2) of m 1 (/i 1 ). 



0.2 




Fig. 2: Monte Carlo profiles for K/K c = 0.001 and h x = 0.01 for lattices of size 128x512, 
256x1024, 512x2048, and 1024x4096 (solid lines from below to above) compared with the nu- 
merical evaluation of the exact result (|B4|). The latter holds for the semi-infinite system. It is 
clearly visible that the Monte Carlo data approach the exact result for increasibng system size. 
The dotted line represents the asymptotic (z — > 0) behavior expressed in Eq. (JsJ) ~ z In z that 
becomes wrong for large distances. 



C. Monte Carlo results 



First we discuss a set of profiles which were obtained with N = 512 by setting h\ = 0.01 
and K/K c varying between 0.996 to 1.004, in steps of 0.001. The data are depicted in 
Fig. 3. Depending on the temperature, we averaged over 10 000 to 30 000 configurations. 
Especially the shape of the critical profile, marked in Fig. 3, is consistent with the scaling 
analysis of Sec. |II B| . It increases up to z ~ 60, then has a maximum and farther away from 
the surface it decays. With increasing distance the influence of the second surface (here at 
z = 2047) becomes stronger, such that the profile has a minimum about halfway between the 
boundaries. (The data of Fig. 3 were not symmetrized after the average over configurations 
was taken.) 

For T above T c (curves below the critical profile), the maximum moves towards the sur- 
face and the regime of growing magnetization becomes smaller. This is consistent with the 
scaling analysis of Sec. [II A] , in this case the growth is limited by the correlation length £. On 
the other hand, for T < T c , the tendency to decay in between the surfaces becomes weaker, 
as the the bulk value of the magnetization in the ordered phase grows. As said in Seep A 



the short-distance growth below T c is described by z Xl x *, where the difference X\ — xs takes 



also the exact value 3/8 f24|| . This time, there is no logarithm however, and the growth is 
steeper than above T c . 
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Fig. 3: Order-parameter profiles for various temperatures below and above T c for fixed h\ = 0.01 
and N = 512 compared with the critical profile. The three curves above the critical pro- 
file were obtained with K/K c = 1.001, 1.002, 1.003. The four profiles below correspond to 
K/K c = 0.999 . . . 0.996. The data are not symmetrized. 

Fig. 4 shows the MC results for the critical point, for different values of h\ in double- 
linear representation. In Fig. 5 the same data are plotted double-logarithmically. The lower 
dashed line shows the short-distance behavior ~ z 3 / 8 In z according to Eq. (|j) and (as already 
discussed in connection with Fig. 2) the MC profiles confirm the scaling analysis of Sec. |llB| . 
The upper dashed line is the pure power law z _1//8 that describes the decay in the regime 
where z is larger than li but still smaller than L or £. In our simulation this regime is only 
reached for relatively large values of hi > 0.5. The uppermost profile (hi = 1.0) obviously 
goes through this regime for 10 < z ^ 60. The finite-size exponential behavior (see Sec. |II C|) 
can be observed in all curves for z > 100. The data of Figs. 4 and 5 are symmetrized, i.e., 
after averaging over configurations we computed the mean value of the left and right halfs 
of the system. Hence, the profiles are only displayed up to halfway between the boundaries 
(here z = 1023). 

The location of the maximum of m(z), z max , as a function of hi is depicted in Fig. 6. 
The maximum z max was determined from the profiles by a graphical method. Error bars 
are estimated. For small values of hi the near-surface growth is limited by finite-size effects. 
Up to about hi = 0.03, the value of 2 max is roughly independent of h\. For larger values 
of hi, z max moves towards the surface. As indicated by the dashed line, the dependence of 
2 max on hi is completely consistent with li ~ h{ 2 obtained from the scaling analysis (see 
(|5bD). For hi = 1.0 (upper curve) the maximum value of m(z) is at the boundary, and the 
magnetization monotonously decays for z > 0. 
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Fig. 4: Monte Carlo profiles for N = 512 at T = T c for hi = 0.005, 0.01, 0.02, 0.03, 0.07, 0.1, and 
1.0 (from bottom to top). The data are symmetrized and are only shown up to z = 1023. 
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Fig. 5: The same profiles as in Fig. 4 in double- logarithmic representation, pronouncing the short- 
distance behavior. For small hi (lower curves) the growth of m(z) described by Eq. @ is clearly 
visible For h\ = 1.0 our Monte Carlo result (upper solid line) is in accord with 771(2) ~ z -1 / 8 
(dashed line). The behavior of the Monte Carlo data for z > 100 is described by the (finite-size) 
exponential decay and corrections to the semi-infinite profiles due to the second boundary become 
stronger. 
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Fig. 6: The location of the maximum z max of m(z) in dependence of h± as obtained from the 
data of Fig. 4 and other profiles for the same system (not displayed). The dashed line represents 



h; 



For larger values of hi we have z D 



l\. For small h\ where l\ becomes larger than 



the finite-size scale L, z max is determined by L. 




Z 

Fig. 7: Order-parameter profiles for K = K c and h\ = 0.01 for system size N = 128, 256, 512, 
1024, and 2048. 



In Fig. 7 profiles for fixed h± and T = T c for various system sizes (ranging between 
N = 128 and 2048) are displayed. With increasing N, the maximum keeps moving away 
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from the surface, z max ~ in the largest system. This means that with these parameters 
we are still in the regime where L < / x . When hi = 0.03 is taken instead, the situation is 
different. The respective MC data are depicted in Fig. 8. For the small system we have again 
increasing 2 max . For the two largest systems iV = 1024 and 2048, however, the maximum is 
roughly located at the same distance from the surface, signaling that now l\ < L. 




0.05 - 



1 10 100 1000 10000 

z 

Fig. 8: Profiles for hi = 0.03 and otherwise the same parameters as in Fig. 7. 



IV. DISCUSSION 

We studied the short-distance behavior of order-parameter profiles in the two-dimensional 
semi-infinite Ising system at or above the bulk critical temperature. Our main goals were a 
detailed understanding of the near-surface behavior of the order parameter and a complete 
scenario for the crossover between "ordinary" and "normal" transition. 

With regard to the short- distance behavior, especially for small surface fields hi, the 
Ising model in d — 2 turned out to be quite special. The functional form is quantitatively 
not captured by the analysis of Ref . || , where the main emphasis was put on the situation 
in d = 3. Here we demonstrated by means of a scaling analysis that a small hi induces a 
magnetization in a surface-near regime that grows as z 3 ^ 8 In z as the distance z increases 
away from the surface (see Eq. (||)). Our findings are not only consistent with available 
exact results (to which we made a detailed comparison), but they also allow a detailed 
understanding of the physical reasons for the growth and its quantitative form. The surface 
magnetization mi generates in the region that is (on macroscopic scales) close to the surface 
and that is much more suceptible than the surface itself, a magnetization m(z) much larger 
than mi. The exponent 3/8 that governs the power-law part of the growth is the difference 
between the scaling dimensions yi (of hi) and (of the bulk magnetization). Eventually, 
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as was demonstrated Sec. |IIB| , the logarithmic factor can be traced back to the logarithmic 
dependence of mi on hi. 

The scenario for the crossover between "ordinary" and "normal" developed in Ref. pfl 
and generalized in this work to include the 2-d Ising model is the following: At T = T c a 
small hi causes the increasing near-surface behavior described above. The magnetization 
grows up to distances li ~ K{ 2 and then the crossover to the power-law decay ~ z -1 / 8 
takes place. With increasing hi the surface-near regime becomes smaller, and eventually for 
hi — > oo the length scale li goes to zero, such that the region with increasing magnetization 
vanishes completely and the situation of the "normal" transition is reached. 

For T slightly above T c our scenario essentially remains valid as long as £ > Only 
when z is of the order of £ the crossover to the exponential decay takes place. For £ < li, 
on the other hand, the growth is limited to the region z < £. 

Concerning three-dimensional systems several experiments were pointed out in Ref. ||, 
whose results are possibly related related to the anomalous short-distance behavior P3| , P5| . 
It would be an interesting question for the future, whether similar, surface- sensitive experi- 
ments in two-dimensional systems are feasible. 
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comments. This work was supported in part by the Deutsche Forschungsgemeinschaft through 
Sonderforschungsbereich 237. 



APPENDIX A: EXACT RESULTS FOR THE SURFACE MAGNETIZATION 



The boundary magnetization mi(hi) was calculated exactly by Au-Yang and Fisher [18 
for the semi-infinite 2-d Ising model at the critical point, in the critical regime, and for a 
strip of finite width at the bulk critical point. Whereas above two dimensions a free surface 
near criticality is paramagnetic in the sense that the response of mi on hi is linear, it was 
found that in d = 2 the function m 1 (/i 1 ) is modified by logarithmic terms. 

In the 2-d Ising model the surface magnetic field hi enters the analysis only through the 
variable 



u 



tanh hi 



which will also play an important role in Appendix [B|. For r 
u « hi) the asymptotic behavior of mi takes the form [IS| 



(Al) 



and small u fin which case 



mi (hi) ~ Bi u {b\ — 2 In u) 



(A2) 



where Bi = 1.5369... and hi = —1.3201... are constants. There is an exact formula 
available for mi (hi) (see Eqs. (2.10) and (2.11) of Ref. |y|) but the asymptotic result flAll) 
was sufficient for carrying out the scaling analysis in Sec. U~E . In Sec. |111 B , we compare our 
MC data for mi with the exact result. The lesson learned from flA~2j ) is that in d = 2 and in 
semi-infinite geometry the dependence of m 1 and hi is not simply linear but modified by a 
logarithmic term: For small hi the surface magnetization behaves as hi In hi. This was not 
taken into account in Ref. HI, where we focussed on the short-distance behavior in d = 3. 
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In the critical region, for r > 0, the leading contribution to the surface magnetization is 
given by 1 18] 



m 1 (h 1 , t)=B iU [ln(2A c r) + B{t/u 2 )\ . (A3) 

This tells us that slightly above the critical point the dependence of m\ on hi is indeed linear, 
but, different from the situation in d > 2, there is a (logarithmic) dependence on the tem- 
perature t. The scaling function B approaches a constant for t/u 2 —>■ and asymptotically 
goes as ~ \n(r/u 2 ) for t/u 2 — > oo, leading back to (|A2| ) . 

In Ref. [|ll| the function mi(hi) was also calculated for a strip of finite width iV and 
infinite lateral extension at the (bulk) critical point. This result is of interest for the com- 
parison with our MC data. Similarly to the case r > 0, the logarithmic term in the limit 
u — > is absent in the strip. The leading contribution can be written in the form 

mi (hx) = Bmibl + cl A" 1 + In A), (A4) 

where b\ = 1.0731 ... and = y/2/2. Thus, as in ( |A3| ) the dependence on hi is linear, with 
an A-dependent amplitude ~ In A. Taking this into account in the scaling analysis, it leads 
to an A-dependent amplitude in the profile, an effect that is clearly visible in Figs. 2, 7 and 
8, where the profiles for different system sizes were plotted. 



APPENDIX B: EXACT RESULTS FOR THE MAGNETIZATION PROFILE 



In a work by Bariev ||19|| , the local magnetization of the semi-infinite 2-d Ising model 
with surface magnetic field was represented in terms of a series of multiple integrals. For 



the system we consider p7[ this solution can be written in the form 



m(n, A, hi) = (l - S 2 ) F{2n/t, l x /£) , (Bl) 

where 5* = sinh(2A), the integer n counts the distance (number of spin columns) from the 
boundary starting with n — 1, and £ is the correlation length defined by 



21-5 



(B2) 



The dimensionaless variables K and hi were defined in @. 

The length scale l x occuring in ( |B1| ) represents the distance up to which the short-distance 
behavior (B) is found. It is explicitely given by |19| 



tanh(A) 1 



(B3) 



The main result of Ref. [0] is the scaling function J-"(x, y) in ( pi| ) expressed in form of a 
series of multiple integrals, with distinct solutions above and below the critical temperature. 
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Since we are interested in the case T > T c , we here only show the result for the former case. 
The solution for the scaling function takes the form [[HJ 



y 



exp(-x/2) exp ^ - f k (x, y) 



(B4a) 



\k=l 



with 



fk(x,y) 



-r / dwi ... dw k 

i K Jo Jo 

ch(wi) 



n ex p 

1=1 

1 ch(wi) 



xch(wi)] x 



ch(wi) + ch(wi + i) ch(wi) + 1 + y 1 



w k+1 = Wi 



(B4b) 



where we abbreviated ch = cosh. 

Several explicit asymptotic solutions were given in Ref. The case z < l\ < £ was 

not considered, however. Also we ourselfs were not able to obtain an explicit solution in 
the limit where £ goes to infinity by zjl\ remain at some finite value, the case, where the 
anomalous short-distance behavior discussed in Sec. 11 B shows up most clearly. The so far 
only available explicit exact solution that is consistent with our findings is given by Eq. (28) 
of Ref. 



0| . It says that, asymptotically for x <C 
function (|B4a| ) is given by 



1 and in the limit y — > + , the scaling 



•F(x,0+) = E 3 | In 0/4) + 1| x 



3/8 



l + -x + 0(x 2 \nx) 
4 



(B5) 



with the constant E% = 0.6987 . . . 

The asymptotic result ( |B5| ) is in agreement with our scaling analysis, which itself provides 
a deeper understanding of the functional form of (|B5| ). Firstly, the exponent 3/8 occurs and 
can be attributed to the difference of scaling dimension as discussed in Sec. [11 B| . Secondly, a 
logarithmic factor shows up, now with an argument ~ z/£ instead of ~ z/h, however. That 
effectively £ replaces l\ for l\ — *• oo can be understood on the basis of our scaling analysis 
and the results reported in Appendix|A], since, as Eq. ( |X3| ) tells us, there is no logarithmic 
dependence of mi on hi for hi — > and finite £. 

In order to obtain profiles for the whole crossover region between the fixed points, one 
has to compute the multiple integrals occuring in ( |B4b| ) numerically to some finite order in 
the series. The series converges rather quickly such that in most cases the termination after 
the third term is sufficient. Only very close to the surface corrections due to higher orders 
become appreciable. 

We evaluated ( |B4D numerically for K = 0.999 K c . The corresponding correlation length 
is £ = 567 (lattice spacings). In Fig. 9 the results for the order-parameter profiles with the 
same values of hi as used in Figs. 4 and 5 are displayed in double-logarithmic representation 
(z = n — 1). Again the power law ~ z^ 1 ^ 8 is displayed for comparison. The resemblance 
between the MC results shown in Fig. 5 and the curves of Fig. 9 is striking. One has to bear 
in mind, however, that in the former case the crossover to the exponential decay was caused 
by finite-size effects. In Fig. 9, on the other hand, it is attributable to the finite correlation 
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length. 



1 

m (z) 



0.1 




1 10 100 

z 

Fig. 9: Order-parameter profiles calculated from Bariev's solution Eq. ( B4| ) for the semi-infinite 
2-d Ising model, with K/K c = 0.999, and for the same values of h\ as in Fig. 4. 
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